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ABSTRACT 

In situ measurements of ion charge states can provide unique insight into the 
heating and evolution of coronal mass ejections when tested against realistic non- 
equilibrium ionization modeling. In this work we investigate the representation 
of the CME magnetic field as an expanding spheromak configuration, where 
the plasma heating is prescribed by the choice of anomalous resistivity and the 
spheromak dynamics. We chose as a test case, the 19 May 2007 CME observed 
by STEREO and ACE. The spheromak is an appealing physical model, because 
the location and degree of heating is fixed by the choice of anomalous resistivity 
and the spheromak expansion rate which we constrain with observations. This 
model can provide the heating required between I.IRq and earth orbit to produce 
charge states observed in the CME flux rope. 

However this source of heating in the spheromak alone has difficulty account- 
ing for the rapid heating to Fe^"^^"'" at lower heights, as observed in STEREO 
EUVI due to the rapid radiative cooling that occurs at the high densities involved. 
Episodes of heating and cooling clearly unrelated to spheromak expansion are ob- 
served prior to the eruption, and presumably still play a role during the eruption 
itself. Spheromak heating is also not capable of reproducing the high Fe charge 
states (Fe^^"*" and higher) seen in situ exterior to the flux rope in this CME. Thus 
while the spheromak configuration may be a valid model for the magnetic topol- 
ogy, other means of energization are still required to provide much of the rapid 
heating observed. 



1. Introduction 

The eruption mechanism of coronal mass ejections (CMEs) is currently an extremely 
active area of research. A variety of mechanisms have been discussed in the literature. 
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Most of these require some form of magnetic reconnection during and after the eruption. 
The associated thermal heating should be an important diagnostic of the eruptive process, 
including the eruption dynamics and ejecta geometry, but until now has been relatively 
unexploited. Such heating gives rise to UV- X-ray emissions during the eruption, that can 
be remotely sensed, and also leaves its imprint in the charge states of ions detected in situ 
near 1 AU, as their ionization state responds to the newly heated plasma. These techniques 
give complementary insights into the heating and evolution. Spectroscopy is in some sense 
more direct, being a prompt signature of the eruption, but with existing instrumentation 
has only been possible relatively close to the solar surface where the signal is strong. Ion 
charge states, on the other hand, are routinely collected and CMEs, and their constituent 
parts (forward shock, flux rope, etc) can be easily identified. Further interpretation of such 
data require modeling of the non-equilibrium ionization (NEI) balance as the CME plasma 
expands outward into the ambient solar wind. In this paper we explore the consequences of 
a different source of thermal heating, that due to anomalous resistivity within a sphcromak 
solution for the CME geometry as it expands into the solar wind. We compare the predictions 
of such a model to STEREO and ACE observations of the 19 May 2007 CME. The structure 
of the paper is as follows. We describe the status of our understanding of ion charge states 
observed in CMEs in more detail in section 2. Section 3 describes the 2007 May 19 CME, 
and some of our motivation for exploring a spheromak model. Section 4 gives the theory 
of spheromak resistive expansion, with some discussion of just how closely the spheromak 
model used here may represent a real CME. Section 5 outlines the NEI model incorporating 
the spheromak resistive heating, and we conclude in section 6 with a summary of our results 
in the context of previous research and discuss future modeling improvements. 



2. NEI modeling of CME charge states 

Interplanetary CMEs or ICMEs have many well-known signatures. Of these, in situ 
measurements of ion charge state distributions hold a unique potential for diagnosing the 
conditions throughout the eruption. The charge states evolve through ionizations and re- 
combinations as a function of temperature and density up to heights of 3 to 6 solar radii 
heliocentric distance but are "frozen in" thereafter. Since low-Z ions "freeze-in" at higher 
densities (i.e. smaller heights or radii) than high-Z ions, they provide complimentary con- 
straints on the temperature evolution. 

A number of studies have examined the existing solar wind composition data and made 
inferences on freeze-in temperatures based on computed ionization distributions appropriate 
to coronal equihbrium (e.g. Zurbuchen and Richardson 2006, Zurbuchen 2004, Lepri and 
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Zurbuchen 200 4, Lepri et al. 2001, Henke et al. 2001, Gloeckler et al. 1998, Henke et al. 
1998). However iRakowski. Laming fc Lepri I ( 120071 ) were the first to conduct time- dependent 
modeling of tlie ion cfiarge state distributions of various elements detected in situ in order 
to draw quantitative conclusions r egarding the thermal input a r id ini tial conditions in the 
corona during the CME eruption. iRakowski. Laming fc Lepri I (j2007l ) modeled the charge 



states of 8 ICMEs from the event list of 
velocities. 



Lynch et al.l (120031 ). chosen to cover a range of 



Rakowski. Laming fc Lepri I ( l2007l ) considered a simplified geometry for the CME ejecta 
and parameterized the evolution. The model velocity and expansion were based on the ob- 
served phenomenology from plane-of-the-sky events showing that typically CM E height-time 



evolution consists of three phases: initiation, acceleration and propagation (IZhang et al. 



200ll . 120041 : IZhang fc Derell2006l : ISui fc Holmanll2003l : iLin et al.ll2005l : ISui et al.ll2005r ). Heat- 
ing of the CME core plasma during the acceleration phase was explored assuming a heating 
rate for the CME plasma proportional to the rate of kinetic energy increase, i.e. a constant 
fraction QE/KE during the acceleration up to a final velocity Vf. In 5 out of 8 sample CMEs 
studied, the dominant Fe charge states were neon-like (16+) or higher, indicating that high 
temperatures, comparable to flares (~ lO^K), are involved. Starting the plasma from this 
temperature and allowing the ions to recombine as they expand could often account for the 
Fe ionization balance, with peaks around Fe^^"*" and Fe^+ (the Ne-like and Ar-like charge 
states, which have small recombination rates to the next charge states down, hence popu- 
lation "bottlenecks" here). However the lower-Z elements placed a limit on the maximum 
starting temperature (at least if assuming ionization equilibrium in the seed plasma). Above 
~ 2.5 X lO^K, O would have been mainly O^^ instead of O^"^ and O^^ as observed and 
would not recombine significantly during the CME evolution. Evidently plasma must start 
out much cooler, and be further heated as the CME accelerates. Thermal energy inputs on 
the order of 2 to 10 times the kinetic energy were needed in 7 out of 8 events studied. 



Similar conclusions to iRakowski. Laming &: Lepri I ( 120071 ) about the thermal energy in- 
put to CMEs have been reached from analysis of ultra-violet spectra taken by SOHO/UVCS. 
Akmal et al.l (l200l[ ) studied a 480 km s"^ CME observed on 1999 April 23, and find a th ermal 
energy comparable to the bulk kinetic energy of the plasma. ICiaravella et al.l (l200l[) give 



similar results for the 260 km s ^ 1997 December 12 CME. More dramatically, iLee et al. 



( 12 0091 ) studied the 2001 December 13 event, and using the combination of [O V] density di 
agnostics and non-equilibrium ionization modeling of O VI fou nd that 75% of the magnetic 



energy must go into heat to match the UVCS observations. iRakowski. Laming fc Lepri 



(120071 ) found a heat to kinetic energy ratio of at least 6 for the same event. 



The motivations behind the current work are twofold. First, to test if our non-equilibrium 



-4- 



ionization modeling of the charge state distributions can hold up to the scrutiny of additional 
constraints on the eruption dynamics coming from STEREO observations in the 2007 May 
19 CME, as well as observations regarding the charge states closer to the sun. Secondly we 
explore whether a more physically and theoretically specified heating model when incorpo- 
rated into our non-equilibrium modeling can reproduce the charge states of a real event with 
reasonable physical parameters. 



3. The 2007 May 19 CME 



The 2007 May 19 CME studied in this paper is qualitatively different from those CMEs 
discussed above. As one of the first major CME eruptions in the STEREO era the 19 



May 2007 event ha s been intensely studied by mu l tiple authors (IGopa. 



Veronie et al. 


I2OO8 




Li et al. 


2006 


]\ 


Liu et al. 


2008 


Bone et al. 


2009: Kerdraon et al. 


2OI0I). It began 



20081 : iKilpua et all 12009 



■swamv et a l. "2OO9I; 



Liewer et al.. ,2009.: 



region AR10956 and was detected as an ICME by STEREO B, ACE and possibly STEREO 
A on 2007 May 22. Multiple heating and cooling episodes were seen in the two days prior 



to eruption (ILiewer et al 



20091 : iBone et all l2009l : iLi et al.ll2008r). Th e total unsigned flux 



declined 17% in the two days prior to the first eruption (ILi et al.ll2008l ). During which time 
at least 4 heating (and cooling) events happened which heated the filament to 1 MK or 
higher ( IBone et al.ll2009[ ). According to lBone et al.l (|2009[ ). the formation and merger of the 
filament which eventually erupted represents 6% of the total unsigned flux change during 
this two day period. 

Our model of the heating during the CME eruption and evolution through interplanetary 
space will need to match both the charge states seen in the filament during eruption and 



Liewer et al. 



detect ed in situ, as well as the density and speed of the ejected material 

J2009h present the STEREO EUV images showing the heating of the filament in the hours 
and minutes immediately preceding the eruption. Hot spots are seen in the EUVI 17lA 
filter (mainly Fe IX and XII) at heights near 1.07 Rq. The magnetic cloud is determined 
to have mainly interacted with STEREO B, which penetrated the center during most of 
May 22 (DoY 142) , while STEREO A only passed through the magnetic cloud periphery 
(IKilpua et al.l 120091 ). The in situ Fe charge states seen with PLASTIC on STEREO B and 
A are shown in Figure 1 0. Charge state data from ACE for multiple ions is shown in Figure 
2. The highest Fe charge states, Fe^^"*", Fe^^"*", are only seen in ACE and STEREO A which 
were not centrally located in the event, but rather on its flanks. Furthermore the appearance 



^ Links to the level 2 data are available at http://fiji.sr.unli.edu/ 
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of the highest charge states is after the cloud passage as defined in iLiu et ahl (|2008[ ). 



The shock front around the 2007 May 19 CME launched at 958 km s~^, the average 
speed in transit to 1 AU was around 700 km s~^ and the sp eed of the MC as it passed 
the STEREO B spacecraft was 482 km s~^ ( iKilpua et al.ll2009l ). The filament that erupted 
from the active reg ion at about the same time only rose at an average speed of 103 km s~^ 
( jLiewer et al.l 120091 ) . but it could have been dragged out by the CME and the fast solar 
wind to reach a similar coasting velocity as the ICME. It had no distinct initiation or 
ac celeration phase, and clearly does not fit within the phenomenology described above and 



m 



Rakowski. Laming &: Lepri I ( 120071 ). which leads us to consider a spheromak model for the 



CME evolution. 



4. Resistive Spheromak model of CMEs 
4.1. Spheromak Field Configurations 



Spheromaks are well known force-free configurations of plasma satisfying condition 
V X B = aB with spatially constant a. They are solutions of the Grad-Shafranov equa- 
tion in spherical coordinates wit h the poloidal current being a linear function of the flux 
(jChandrasekhar &: Kendall! 119571 ). The basic spheromak solution is 



B, = 25o^^i-^cos^ 
ar 

5, = -V-iM±^!lM3in^ 

ar 

B^ = Boji (ar) sin 9 



where the spherical Bessel function ji can be expressed in terms of elementary functions, 
j'i(x) = sinx/x"^ — cosx/x. The parameter a is related to the size of the spheromak R, 
defined by the surface where radial magnetic field is zero, given by the solution of ji = 0; 
R = Ca/a with = 4.49. 

The spheromak model of CME topology is an attractive alternative to the usual fiux 
rope model. It is self-contained with no unknown "length of fiux rope" parameter to main- 
tain connectivity to the sun and the heating within it is mathematically specified once the 
anomolous resistivity is chosen . The torus structur e, seen on ed ge, may even resemble the 
f ilame nt shape in some CMEs (jKataoka et al.ll2009l ). Previously, iLyutikov &: Gourgouliatos 
( I2OIOI ) found ideal self-similar solutions for expanding spheromak, with electric fields 



r a 

E = — e^ X B 

c a 



(2) 
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(dot denotes differentiation witli respect to time) and tlie corresponding non-radial velocity 



In tliis Section we discuss continuous lieating of an expanding splieromak, generalizing tlie 
previous analysis to non-ideal self-similar expansion with E • B 7^ 0. 

Physically, the parallel component of the electric field should be related to the cur- 
rent density through Ohm's law. Formally, the procedure described below breaks down the 
assumption of self-similarity, since the value of the parallel electric field component is not 
linearly proportional to the current density. Still, we assume that resistivity plays a sub- 
dominant role, so that the expansion remains approximately self-similar. In other words, we 
assume that resistivity leads to small deviation from the ideal self-similar expansion. The 
resulting dynamics remain self-similar, by assumption, but are a little different from the 
ideal case. We are interested not in the detailed properties of local resistive heating, but 
in general scaling relations. As we show below, the requirement that an expanding CME 
dissipates some of its initial energy in the form of heat, detectable by its effect on element 
charge states by the time it reaches Earth orbit, can be used to estimate the anomalous 
resistivity. 

Expressing electric field in radiation gauge E = —dtA/c, where A is vector potential, 
we find from the force-free condition. 



This equations highlights two important points. First, it shows how the time-dependences 
of the poloidal components of the vector potential are related to the toroidal components, 
Ag,Ar oc a{t)A^ (this is related to the conservation of poloidal and toroidal fluxes). Second, 
we get one equation for the two functions Ar,Ag. Since spheromak solutions are linear, a 
general solution is a linear combination of the two solutions. The two particular solutions, 
which in a static case correspond to two equivalent choices of the poloidal components of 
vector potential, in a time varying case correspond to different physical processes. Accord- 
ingly, there are two types of solutions, corresponding to two choices of vector potential in 
Eq. dl]), 7^ and Ag ^ 0. We consider them in turn. 



field 




(3) 




(4) 
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4.2. Resistive expansion of first type 



Let us assume that scale with time as an arbitrary function of the expansion param- 
eter, B(f, oc /(a/ao). The first type of solutions has Ar = and, from equation 1, 



, , , (sin(ar) - r —dz) 
Ae = -Bof{a/ao) ^ ^ ' ^ ^sin^. 



(5) 



The first electromagnetic invariant E • B oc [af — 2a f^. Thus, if / = {a/aoY the flow is 
ideal, consistent with the scaling chosen in lLyutikov fc GourgouliatosI (l2010[ ). More generally, 
if we allow / = (a/ao)^'*'"^, we find the parallel electric field 



El, = m'BUa/aor"'^^^F,ira,i 



(6) 



where Fi{ra, 9) is a lengthy function not given explicitly here. The case of m = corresponds 
to ideal expansion. Since we expect that reconnection eliminates magnetic flux, physically 
realizable solutions correspond to m > 0. Magnetic flux decreases as a™ oc R~"^. 



4.3. Resistive expansion of second type 

The second type of dissipative expansion corresponds to = and 

Ar = Borf{a/ ao) ji cos 9 (7) 
The parallel electric field in this case 

El = B',ia/aor^"'~'^^F,{ra,9). (8) 

For this solution the resulting electric field has a component along magnetic field for any 
m. The resistive effects in these solutions are induced exclusively by spheromak expansion. 
Note, that this is different from resistive decay of a stationary spheromak, which proceeds 
homologously. Thus, expansion of the spheromak generally leads to appearances of parallel 
electric fields and associated dissipation, which can proceed much faster than resistive decay 
of stationary spheromaks. 



4.4. Expansion-driven dissipation 



Let us concentrate on the resistive expansion of first type, which produces equatorial 
dissipation concentrated in the flux rope, as in the CME observations we attempt to model. 
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The dynamical models described above relate the resistivity to the overall djTiamics. The 
volume integral over En (see Eq. (Q) is 



„,2 / D \ 2m+l 

Ef^dV = 0.033m'B',Rl^ ( ^ ) (9) 



which leads to the dissipation rate 

-Adv = omSmB^Rl ( ^ I t; (10) 
r] \ R J R 

with the identification of the resistivity rj = mRv/c^. This is justified by integrating the 
energy dissipation rate above between _Ro and R to give a total dissipated energy 



""Ur = omS^^vRt (l-( ^V™ 1 OmS^^vR' X 2mln ( ^ ] (11^ 



Ro 



V 



where the last step assumes {R — Rq) / Rq « 1- With total spheromak (magnetic) energy 
[/ = 1.57 X 10-^ B^R'^/R X {Ro/Rf"", the difference in energy between Ro and R is SU = 
U (-Ro) (l ~ (-RoZ-R)^™^^) ; of which a fraction 2m/ (2m + 1) is dissipated. Writing 

6U ^ 2mU (Ro) In ( ^ ) = OmSB^Rl In ( ^ ) (12) 



and equating with ^dr in equation 11 yields rj = mRov/c^. This is derived assuming an 
infinitesimal expansion. The generalization to arbitrary expansions gives 

f] ~ mRv/(?. (13) 



In the second case the dissipation rate 

// p \ 2(m+l) 
¥.-UV = BlRl{-^\ VF.im) (14) 

where F2{m) ^ 0.033m — 5 x lO^'*, very similar to the first case when integrated over the 
spheromak volume. The distribution of the heating is different in the two illustrated 
in Figure [31 The amount of energy available for heating therefore depends on the value of 
the anomalous resistivity and varies with position in the spheromak about an average value. 
The distribution of volumetric heating rates for the first type of dissipative solutions (of 
most interest here; see below) is given in Figure ID 

To what extent can a spheromak approximate a real CME? Flux rope structures, similar 
to the equatorial portion of a spheromak are frequently observed in situ, through the rotation 
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of the magnetic field as the spacecraft penetrates the CME plasma. The precise origin of the 
fiux rope is not clear. Many authors assunie that such a structure eme rges fully formed from 
the photospher e fe.g..lFan &: GibsonI 20041: iMagara &: Longcopdl2003l ). Altern atively in the 



break out, (e.g. lAntiochos et al. 



2OOII ) or fiux cancelation (e.g. [Linker et al. 



1999 



Lynch et al.ll2004[ ). tether-cutting (e.g. iMoore et al. 



2003! ) models the fiux rope forms as the CME 



erupts by the reconnection along an arcade of magnetic loops. In such a case, the ends of 
the fiux rope remain attached to the Sun during the eruption. A spheromak however is a 
self contained magnetic field configuration supported by internal currents, with no external 
attachment, and presumably if appropriate, must emerge fully formed from beneath the solar 
photosphere. 

In general, the models in which a fiux rope forms by reconnection predict an acceleration 
phase at the onset of fiare reconnection. By contrast, the 2007 May 19 CME appears 
to have moved out from very low down at almost constant speed. There is no obvious 
feature in its height-time plot that would correspond to the onset of an acceleration phase 
coupled with reconnection to form a fiux rope, suggesting that the observed fiux rope must 
have formed much lower down. Thus its observed trajectory suggests that a spheromak 
might be a more appro priate description here, than it wo u ld have been for other CMEs 
considered previously in iRakowski. Laming fc Lepri I ( 120071 ) . iNakagawa fc Matsuokal (120101 ) 
find a similar trajectory for the magnetic cloud observed by ACE and Nozomi on 1999 April 
16-18. They also conclude that the magnetic cloud is better fitted by a toroidal, as opposed 
to cylindrical, fiux rope, though the model fitted is still approximate. Both fiux ropes and 
spheromaks conserve helicity as they expand, and to d o so must decrease th eir magnetic 
energy, a part of which can end up dissipated as heat. iKumar fc RustI (119961 ) consider an 
expanding fiux rope, and find temperatures reaching 1 — 2 x lO^K. We apply similar ideas, 
derived within the context of a spheromak, to the 2007 May 19 CME. We assume that the 
expansion velocity of the spheromak entering the expression for the heating is given by the 
observed motion of the CME. 



5. Spheromak non-equilibrium ionization model of the 19 May 2007 CME 



As in iRakowski. Laming fc Lepri I ( l2007l ) we modeled the charge states within the CME 
ejecta for a variety of ions using an adapta tion of the BLASPHEMER (BLASt Propaga- 
tion i n Highly EMitting EnviRonment) code (ILaming fc Grunll2002l l2003l : iLaming fc Hwang 



20031 ) . which follows the time dependent ionization balance and temperatures of a Lagrangian 
plasma parcel as it expands in the solar wind. The fundamental equations are outlined 
below. Further details of the ionization and recombination calculation can be found in 
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Rakowski. Laming &: Lepri I (120071 ) and references therein. 

The density riiq of ions of element i with charge q is given by 



Hq 



dt 



(C*io?i,(j— l^i q~l CiQ^^q^iq^ ~\~ fie ~l~ Cdr,q+l) q+1 {Crr,q ~l~ Cd^ q^ 



n. 



iq 



(15) 

where Cion,q, Crr,q, Cdr,q dxe the rates for electron impact ionization, radiative recombination 
and dielectronic recombination respectively, out of the charge st ate q. These rates a re the 
same as those used in the rec ent ionization balance calculations o f iBryans et al.l (120061 ) . with 
more recent updates given in iRakowski. Laming fc Lepri I (120071 ). The electron density rig is 
determined from the condition that the plasma be electrically neutral. The ion and electron 
temperatures. Tig and are coupled by Coulomb collisions by 



dTi 



iq 



dt 



-0.13n, 



{Tig-Te)q'n,J{q + l) /In A 



MigTe 



3/2 



{J2iq ^^^<l) 



37 



(16) 



and 

dTe 
dt 



O.lSne j2 C^^i - 'l^'^^il (9 + 1) /In A 



3/2 



iq 



iq 



E 



iq ^iq 



37 



Te f dUe 
He V dt 



2 dQ 



+ 



3ne/cB dt l.briek^ 



(17) 

Here Miq is the atomic mass of the ions of element i and charge q in the plasma, and In A ~ 28 
is the Coulomb logarithm. The term in dQ/dT represen ts plasma energy losses due t o 
ionization and radiation. Radiation losses can be taken from [Summers fc McWhirterl ( 1l979l ). 
In a CME from a filament eruption the densities and temperatures are such that radiative 
losses can be important, unlike most applications in the solar wind where they are generally 
negligible. The term — (Te/rig) [drie/dt)-^^ gives the reduction in electron temperature when 
the electron density increases due to ionization. Recombinations, which reduce the electron 
density, do not result in an increase in the electron temperature in low density plasmas, since 
the energy of the recombined electron is radiated away (in either radiative or dielectronic 
recombination), rather than being shared with the other plasma electrons as would be the 
case for three-body recombination in dense plasmas. We include the last electron heating 
term to model the Ohmic dissipation given by equation 10 for the heating in the spheromak. 

In addition to ionization and recombination there is also the geometry to consider in 
the density and temperature evolution. For lack of a better motivated profile we chose a 
uniform initial density and temperature, exploding from an initial radius of [zq — I) where Zo 
is the starting height from the center of the sun in solar radii. Thus the underlying density 
evolution for adiabatic expansion goes as [{z — l)/{zo — 1)]^, which goes over to (z/zq)'^ as 
z » zq and the CME expansion more closely follows that of the ambient solar wind. 
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There are two questions we are trying to answer with our simulations. First, how well 
can a spheromak solution explain the charge states seen in the 19 May 2007 event? To this 
end, two sets of models were explored. In model 1 we attempted to match both the filament 
heating and the ICME charge states using only the heating due to anomalous resistivity in 
the spheromak. In model 2 we allowed the plasma to already have been heated during the 
eruption to temperatures that would explain the EUVI observations, but then followed the 
subsequent heating in the spheromak. Second, given that fast CMEs, particularly at solar 
maximum, often show Fe charge states as high as Fe^^+ in situ, what are the highest Fe 
charge states obtainable by anomalous resistive heating in a spheromak model? In other 
words, how general could such a model be? 

The baseline model 1 starts the plasma at a temperature of lO^K, a density of 10^°cm~^, 
and a velocity of 10 km s~^ accelerating to 500 km s~^. These were chosen as appropriate 
temperatures and densities for cool filament material, and the evolution is consistent with 
the final density seen in situ. We initially explored the effect on the final charge states 
of changing the anomalous resistivity, and thus the amount of magnetic flux dissipation, 
through the parameter m. We assumed heating at the average rate for the spheromak 
volume. Values of m of 0.01 and below leave the charge states essentially unaffected by 
the spheromak heating. (The charge states still evolve with the adiabatic expansion which 
changes the temperature and density as the CME erupts). At m values between 1.5 and 
2 one reaches the limit of diminishing returns whereby increasing m does not result in 
higher charge states, due to the factor (Rq/R)'^'^ in equation 10. This corresponds to an 
anomalous resistivity between 8 x 10~^s and 10~^s, as much as 10^^ times the Spitzer value, 
with a maximum Fe charge state of Fe^'^"'". We note that increasing the final velocity or 
the acceleration rate does not substantially change this limit. A more dramatic spheromak 
expansion starting from a smaller sphere does heat the material faster, but it also reduces the 
density too quickly such that (a) the final densities are substantially lower than seen in situ 
and (b) ionization declines such that the highest charge states are not produced. Thus we 
conclude that the high charge states over Fe^^+ in fast CMEs cannot be reproduced purely 
with the spheromak solution and anomalous resistivity. This leaves the question of how well 
the spheromak can explain this particular 19 May 2007 CME. 

We choose m — 0.75 and compute ion charge states averaged over the whole spheromak 
assuming the distribution of heating rates in Figure 3. Since the 2007 May 22 CME flux 
rope does not exhibit Fe more ionized than Fe^^+, its highest charge states can be explained 
with the basehne model 1 and an average m — 0.75. The overall weighted charge states 
of this model are shown in Figure 5. There are two obvious discrepancies between in situ 
observations and the Fe charge states modeled here. First, there are too many low charge 
states from the unheated portion of the spheromak. This is not a major concern as it is 
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possible that that portion of the spheromak is simply not filled with plasma. Second, Carbon 
charge states seem to indicate that the plasma started hotter than in these simulations. 

Additionally one can consider whether the spheromak heating could explain the EUV 
observations of the filament around launch. Figure 6 shows the early evolution of the charge 
states corresponding to the average heating rate near the sun in the baseline model 1 with 
m = 0.75. This s tarted with the plasma at I.IRq, above the heights of I.OTRq where 



Liewer et al.l ( l2009l ) first see Fe IX and Fe XII emission. Starting simulations below this 
height, with a density that evolves to become the observed density in situ at earth orbit, 
leads to strong radiative cooling that keeps the electron temperature below ~ lO^K until 
a radius of l.li?0 is reached, whereupon the density is sufficiently reduc ed to lengthen the 
radiative cooling time. Hence the heating observed by lLiewer et al.l ( 120091 ). beginning before 
the actual eruption began, must have had a different origin. 

Testing models that started at lower heights, we found that indeed Fe IX and even Fe 
XII could be produced within a short distance of launch. We did not use these for our in situ 
modeling because for reasonable starting filament densities the final density due to the rapid 
expansion was too low. Overall we conclude that while reproducing the EUV observations 
may be possible within a spheromak context, gi ven the sequence o f heating and cooling of 
the filament leading up the event as described by iBone et al.l (l2009l ) , it is also plausible, and 
probably more likely, that whatever mechanism was heating the filament pre-eruption was 
responsible for the filament heating during the eruption. 

The above discussion naturally leads us to our second model to determine if spheromak 
heating starting from temperatures seen in EUV observations can explain the in situ charge 
states. In Model 2 we started the plasma at lO^K, 10^cm~^, and a velocity of 100 km s~^ 
similar to the velocity of the filament from (ILiewer et al.ll2009l ) allowing it to then accelerate 
to the coasting velocity of 500 km s~^. Starting from this high temperature, the additional 
heating of the spheromak has no noticeable effect on the charge states. While the tempera- 
ture was chosen such that the charge states at low heights match observations, the final Fe 
charge states fall short of those detected in situ. We also experimented with starting at lO^K 
but the higher density and/or at a slower velocity. However, radiative recombination domi- 
nated bringing the charge states and temperature rapidly back to their lO^K levels. This has 
implications for the rapidity and amount of energy that had to be dumped into the filaments 
to produce their high charge states. A thorough examination of this phenomenology is war- 
ranted but we leave this to future work. Overall, starting from a higher initial temperature 
in a uniformly filled spheromak did not solve the problem of matching the highest Fe charge 
states. 
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Discussion and Conclusions 



We have explored modeling this filament eruption as due to the heating provided in an 
expanding spheromak. Figure 6 shows the evolution of H, C, O, Si and Fe charge states 
obtained with an anomalous resistivity increased from the classical (i.e. Spitzer) value by 
a factor around 2 x 10^ — 6 x 10^^, as the plasma temperature varies between 10'^ K and 
10^ K. Anomalous resistivity is thought to arise in conditions where plasma turbulence can 
scatter current carrying electrons, and can do so more effectively than the Coulomb collisions 
between electrons and ions. It may be estimated by replacing the isotropization frequency 
by Coulomb collisions with the relevant wave-particle interaction rate. 

Solutions of the one dimensional Vlasov equation desi gned to model the Buneman insta- 
bility in a reconnecting current sheet (IWu &: Huangl 120091 ) . and therefore model electrons in 
Langmuir turbulence, give quiescent values of the anomalous resistivity of xlO^ the Spitzer 
value, with transient values as high as xlO^ in a background plasma temperature of 10'' K. 
In our case, dropping the anomalous resistivity to xlO^° the classical value, the dominant 
Fe charge state becomes Fe^"*". The ratio of the electron plasma frequency to the collision 
frequency at IARq in Figure 6 evaluates to ~ 10^, and proba bly repres e nts th e maximum 
enhancement in resistivity than one may reasonably expect. iLin et al.l (120071 ) estimate a 
much higher resistivity from the obse rved width of cur rent sheets trailing CMEs, as high 
as 10^^ times the classical resistivity. iBemporadI (120081 ) revisits this and argues that what 
is observed as a macroscopic current sheet is in fact an assembly of microscopic current 
sheets, each with wi dth corresponding to anomalous resistivities similar to those given by 
( IWu fc Huangl |2009[ ). although possibly attributed to different modes of turbulence. We 
therefore conclude that magnetic energy dissipated in the expansion of a spheromak may ac- 
count for the filament heating between 10^ K and 10^ K, but is unlikely to be the explanation 
for the high Fe charge states characteristic of temperatures of order 10^ K and higher. The 
nature of the heating specified in the spheromak solution is such that the heating rate due 
to anomalous resistivity cannot be increased without limit, even if an anomalous resistivity 
as high as allowed by the solution could be justified physically. 

Higher charge states, beyond Fe^^"*" must therefore result from other means of heating. 
Noting that in the 2007 May 19 CME, the highest Fe charge states are observed exterior to the 
flux ropes, we speculate that reconnection associated with the eru ption is most likely. T ype 



II and Type HI radio bursts were observed during this eruption (IKerdraon et al 



times when the CME shock was between 1.2 and 1.4 Rq heliocentric distance ( IVeronig et al. 



20101) at 



20081 ). Bearing in mind that the plasma we observe and model will be at lower altitude 
that the CME shock, this corresponds very well to the epochs of heating shown in Figure 
6. Recent results suggest that reconnection might be an efficient means of heating electrons 
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Drake et al. 


2006; 


Oka et al. 


2010) 



In earlier work (IRakowski. Laming fc Lepri 1120071 ). we expressed the heat input to the 
CME in terms of its kinetic energy input by relating the heating to the acceleration, tac- 
itly assuming that the mechanism of acceleration inevitably also supplies heat to the CME 
plasma. In the case of the spheromak, the heating is related to the expansion, but not the 
acceleration, so the ratio of heat to kinetic energy doe s not have a natural interpr e tation . 
We estimate it, though, to facilitate comparison with iRakowski. Laming fc Lepri I (120071 ). 
From equations 11 and 12, 



/ edt OmSB^R^ In (R/Ro) 



\Rq 



which close to the Sun evaluates to ^ 1. Th is is si milar to, but slightly lower than typical 



values found by IRakowski. Laming fc Lepri I (120071 ) . The CMEs studied in that work fre- 
quently had charge states of Fe up to neon-like, rather higher than in the flux rope studied 
here which appear to be consistent with our estimate of thermal to kinetic energy in the 
CME. 



This work has been supported by NASA Contract NNG08EK62I and by basic research 
funds of the Office of Naval Research. 
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Fig. 1. — Fe charge state distributions detected by STEREO A and B during the 2007 May 
19 CME passage. Highest Fe charge states are detected by STEREO A, though the magnetic 
cloud is more clearly detected in STEREO B. Magnetic field and other data indicate cloud 
passage between 3:36 hr and 16:34 hr on day 142 (22 May) for STEREO B and 21 May 19:12 
until 22 May 00:14 for STEREO A Jhiu et al.lboosi ). The dashed vertical white lines show 
the temporal extent of the magnetic cloud (or fiux rope) in this event. 



This preprint was prepared with the AAS E^TJrjX macros v5.2. 
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Fig. 2. — ACE charge state data for C, O, Si and Fe in 8 hr bins. Earliest times have been 
offset upwards for display purposes. The highest Fe charge st ates only appear starting in 
Day 143 after the passage of the interior of the magnetic cloud. iLiu et al.l (l2008[ l report the 
ICME passages through ACE as 21 May 22:19 to 22 May 12:43. 



Fig. 3. — Normalized dissipation rate e as function of radial distance r / R{t) and polar angle 
9 for two self-similar dissipative structures, m=l. Left Panel: solutions ([MH])- Right Panel: 
solutions ([3ED. 



-20- 




0.00 



0.0 0.5 1.0 1.5 2.0 

heating rate/average heating rate 



Fig. 4. — Distribution of volumetric heating rates in the spheromak model of the first type, 
solution ([Mn]), in terms of the average heating rate over the whole volume. 
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Fig. 5. — Final charge state distributions of C, O, Si and Fe for model 1. The individual lines 
show the contributions from regions within the spheromak with different heating rates. The 
histograms show the cumulative charge distribution assuming a volume weighted summation. 
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Fig. 6. — Evolution of charge state distributions for the highest anomalous resistivity con- 
sidered in model 1 for H, C, O, Si, and Fe. Thick lines highhght 0+^ Si+^ Fe+^ Fe+^\ 



